function [E1,E2,E3] = calculate_(h)
True=cos(1);
E1=(sin(1.0+h)-sin(1.0))/h-True;
E2=(sin(1.0)-sin(1.0-h))/h-True;
E3=(sin(1.0+h)-sin(1.0-h))/(2*h)-True;
end
